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摘 要 :针对 压缩 感知 理论 在 宽带 频谱 感知 领域 应 用 时 重 构 精 度 差 的 问题 ,根据 平稳 信号 在 频 域 所 表现 出 的 稀疏 特性 ， 
提出 了 一 种 基于 P-Ifourier(Partial-Inverse fourier) 观 测算 阵 的 宽带 压缩 频谱 感知 方法 。 新 方法 首先 将 频谱 感知 问题 建 模 
为 一 个 典型 的 压缩 感知 问题 ， 利 用 相关 性 能 优良 的 标准 正 交 侍 里 叶 基 构造 观测 矩阵 ， 使 观测 矩阵 具有 良好 的 重 构 性 能 
和 重 构 精 度 。 仿 真 结果 表明 ， 相 比 于 高 斯 随机 观测 矩阵 和 点 入 式 混沌 序列 -循环 Toeplitz 结构 观测 天 阵 ， 该 方法 在 较 低 
信 骂 比 环境 下 能 够 明显 降低 信号 重 构 的 均 方 误差 ， 并且 在 相同 条 件 下 的 重 构 概率 得 到 了 明显 改善 。 
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Broadband compressed sensing method based on P-Ifourier observation matrix 
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Abstract: Aiming at the problem of poor precision in the application of compressed sensing theory in the field of broadband 


spectrum sensing, and according to the sparse characteristics of the stationary signal in the frequency domain, this paper 


developed a broadband compression Spectrum sensing method based on P-Ifourier (Partial-Inverse fourier) observation matrix. 


The new method translated the spectrum sensing problem into a typical compressed sensing problem, and used the standard 
orthogonal Fourier basis observation matrix which has the excellent incoherence performance to build the observation matrix in 
order to have good reconstruction performance and reconstruction precision. The simulation results show that compared with 
the Gaussian random observation matrix and the embedded chaotic sequence - cyclic Toeplitz structure observation matrix, this 
method can significantly reduce the mean square error of signal reconstruction in the lower SNR environment, and under the 
same conditions it can significantly improves the probability of reconstruction. 
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前 ， 许 多 国内 外 学 者 针对 这 一 领域 有 不 同方 面 的 研究 。 
Tian 等 人 首先 提出 并 验证 了 宽带 压缩 频谱 感知 的 有 效 性 。 之 
随 着 信息 时 代 的 不 断 发 展 ， 人 们 对 无 线 电 资源 的 需求 也 日 后 Tian 等 人 中 又 提出 了 一 种 基于 循环 检测 的 宽带 压缩 频谱 感知 
益 增 多 ， 由 于 无 线 频谱 资源 是 有 限 的 ， 因 此 ， 如 何 高 效 合理 算法 。 文 献 [3] 提 出 了 一 种 确定 性 人 藤 入 式 混沌 序列 -循环 Toeplitz 
分 配 和 管理 频谱 资源 是 当前 急需 解决 的 重要 难题 。 基 于 认 知 无 ”结构 的 观测 抢 阵 ， 但 是 该 矩阵 依赖 于 所 选取 的 伪 随 机 序列 的 自 
线 电 理论 的 频谱 感知 技术 通过 不 断 探索 频谱 空 穴 而 动态 地 分 配 。 相关 特性 。 文献 [4] 提 出 了 一 种 优化 的 自 适 应 宽带 压缩 频谱 感知 
频谱 资源 为 解决 这 一 问题 提供 了 思路 ， 但 同样 面临 很 多 挑战 ， 算法 ， 将 观测 矩阵 优化 与 压缩 采样 的 自 适 应 过 程 相 结合 ， 但 是 
挑战 之 一 就 是 随 着 无 线 电信 号 带宽 的 不 断 增加 ， 基 于 传统 奈奈 ”产生 该 类 和 矩阵 构造 过 程 较为 复杂 ， 且 没有 利用 频 域 稀 疏 信号 的 
斯 特 采 样 定律 的 硬件 设备 越 来 越 难以 承受 因 高 采样 率 而 带 来 的 特点 进行 构造 ， 因 此 重 构 效果 有 待 进一步 挖掘 。 文 献 [5] 将 交叉 
载荷 。 针 对 这 一 问题 ， 将 压缩 感知 (CS: Compressed Sensing ) 验证 理论 引入 压缩 感知 领域 中 ， 对 观测 矩阵 进行 优化 ， 但 该 久 
理论 引入 到 宽带 频谱 感知 技术 中 已 成 为 当今 的 研究 热点 之 一 。 车 是 不 断 对 高 斯 随机 秆 阵 进行 优化 得 来 的 ， 因 此 加 大 了 硬件 负 
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所 研究 的 频谱 感知 场景 ,本文 构造 了 P-Ifourier 观测 矩阵 ， 并 将 


民 据 感知 环境 构造 或 
住 恢 复 的 前 提 。 本 文 首先 对 傅 里 叶 算 
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E 阵 具有 理想 的 约束 等 距 特 性 〈RIP: restricted 
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电 和 矩阵 有 较为 明显 的 改善 。 


分 为 工 个 子 频 


互 不 重 登 。 主 用 


户 可 以 随意 地 占 / 


j 这 些 子 频带 ， 


的 子 频带 ， 它 介 
density，PSD) 接近 于 零 ， 医 
[4], 


此 PSD 


] 的 功率 谱 密度 


在 频 域 


设 待 检测 的 宽带 信号 为 xD， 它 可 以 表示 为 


其 中 : 


心 频率 ,n(t) 表 示 


w, 表示 第 i 个子 频带 的 出 


Xx(1) = py QcOs27 ft +n(t) 
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(power spectrum 


上 表现 有 一 定 的 稀 


值 ， f 表示 第 i 个子 频 带 的 中 


可 以 近似 为 零 。 
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性 ， 而 在 频 域 上 呈现 稀疏 特性 ， 
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和 ,在 


定 信 噪 比 条 件 下 ， 
信号 在 时 域 上 都 不 具有 稀疏 特 
此 将 信号 x(t) 表 示 到 频 域 得 
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故 办 ; 正 为 离散 傅 里 时 变换 基 ，s 为 NN 维 信 号 x 在 频 域 下 的 稀 


疏 系 数 向 量 ， 向 量 s 中 只 及 个 系数 不 为 零 (K<N)， 其余 系数 
为 零 或 近似 为 零 ， 称 向 量 s 是 K 阶 稀 琉 的 。 
信号 的 稀疏 表示 是 压缩 感知 理论 应 用 的 第 一 步 。 将 得 到 的 
频 域 稀疏 系数 向 量 * 进行 压缩 采样 可 得 观测 向 量 
y=As (3) 
专访 : y 为 MX1 维 观测 向 量 (M<<N)，4 为 MXN 维 观测 算 
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对 于 压缩 采样 模块 中 的 观测 矩阵 4， 
局 部 傅 里 叶 和 矩阵 以 及 托 普 利 兹 (Toeplitz) 矩 
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稀疏 表示 。 。 稀 琉 系数 向 量 *s ， 压缩 采样 “| 观测 向 量 y 
s=Fx y=As 
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自 相 关 特 性 ， 因 此 应 | 
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E 阵 中 随机 
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于 本 文中 的 频谱 感 和 


傅 里 叶 和 矩阵 只 适 ) 


域 


于 时 


0 场景 ， 因 此 针对 本 文 


在 节 详 细 介绍 该 矩阵 的 构造 思想 及 其 性 能 。 
图 1 及 式 (3) 可 知 ， 只 要 利用 恢复 算法 和 测量 向 量 了 本 
构 出 稀 玻 向 量 s, 就 能 获得 信号 的 频谱 信息 。 对 于 式 (3) 的 求解 ， 
可 以 通过 求解 基于 eo 范 数 最 小 化 的 问题 来 恢复 信号 的 频谱 s。 
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Ss 


min St.y=As (4) 
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此 方法 求 得 的 解 是 最 优 解 , 但 由 于 M<N, 式 (4) 的 求解 属于 
NP-Hard 难题 。 针 对 此 类 问题 ， 压 缩 感 知 理论 中 的 一 个 主要 方 
法 就 是 求解 一 个 基于 el 范 数 最 小 化 问题 并 且 利用 线性 规划 问题 


来 求解 重 构 信号 中。 文献 [10] 证 明了 eo 最 小 范 数 在 一 定 条 件 下 
与 e 最 小 范 数 等 价 。 于 是 式 (4) 可 以 转换 为 
min fl， st. y= As (5) 
要 有 两 类 算法 ， 一 种 是 凸 优化 算法 ， 一 种 


得 到 全 局 最 优 解 ， 但 是 算法 复杂 度 高 ， 速 
迭代 的 思想 , 虽然 恢复 精度 没有 前 者 高 ， 
度 快 等 优点 ， 前 常用 的 一 类 算法 ， 
主要 包括 匹配 追踪 算法 (MP: Matching pursuit)、 正 交 匹 配 追 踪 
算法 (OMP: Orthogonal Matching Pursuit ) 和 压缩 采样 匹配 追踪 
。 其 
前 最 为 常 


三 | 
元 


等 


算法 (CoSaMP:Compressive Sampling Matching Pursuit ) 等 
中 正 交 匹配 追踪 算法 收敛 速度 快 、 恢 复 精度 高 ， 是 
的 恢复 算法 ， 本 文 将 在 仿真 实验 中 运用 此 算法 。 


ey 


2 ”观测 矩阵 的 构造 及 证 明 
2.1 观测 矩阵 构造 条 件 


在 压缩 感知 理论 中 ， 稀 疏 信 号 能 和 否 重 构成 功 取决 于 观测 矩 
阵 是 否 满 足 约束 等 距 特 性 ， 文 献 [11] 对 约束 等 距 特 性 给 予 了 详 
的 介绍 ， 它 可 以 表示 为 
G-G)ls| <|4s), < C+6)|sl, (©) 
其 中 : s 为 任意 阶 稀 琉 向 量 ， 常 数 。 但 是 式 (6) 在 实际 证 
明 过 程 中 很 难 直 接 计 算 ， 因 此 式 〈6) 可 以 同等 转换 为 
1-6<4 A A) SAAA)<1+6, (7) 
其 中 : 7 是 索引 集合 , 满足 Tc{1,2,…NN} 且 |T|<K。AT 是 从 观 
测 矩 阵 4 中 选取 索引 集合 了 所 对 应 的 列 而 构成 的 和 矩阵。 
4 Ca4) 与 1 C4a4) 分 别 表示 和 矩阵 484 的 最 大 、 最 小 特征 
值 ， 由 式 (7) 可 知 : 8. 越 小 ,观测 矩阵 的 最 大 、 最 小 特征 值 越 
接近 于 1， 观 测 矩 阵 的 RIP 性 能 越 优 异 ， 重 构 性 能 越 好 。 
除 此 之 外 ， 文 献 [12~15] 指 出 Donoho 等 人 所 提出 的 相关 性 
判别 理论 能 够 比较 直观 地 判别 观测 矩阵 的 形态 以 及 性 能 。 它 可 
以 简单 的 表示 为 


AS 


L(A)=VN -max | dij | 1<i< J <N (8) 
其 中 : 元 为 观测 矩阵 4 的 原始 方 阵 ，| w ,| 表示 方 阵 了 中 第 i 行 
第 7 列 元 素 的 幅 值 。 简 单 的 说 就 是 原始 方 阵 的 所 有 元 素 中 最 大 


的 幅 值 定义 为 该 矩阵 的 相关 性 ， 风 罗 <[LVN] ,以 帮 值 越 小 ， 
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Vv 合作 
录用 稿 刘 洋 ， : 基于 了 - wb 观 ; 由 YR 作 期 和 和 
观测 矩阵 的 相关 性 越 小 ， 重 构 性 能 越 优 异 。 其 中 : W =e-27*, 算 阵 Fn 的 每 一 列 不 是 单位 向 量 , 因此 Fx 不 
2.2 ”观测 矩阵 的 构造 是 酉 矩阵 。 
针对 频 域 稀疏 的 信号 ， 青 结合 本 文 提出 的 频谱 感知 模型 ， b) 将 矩阵 Fx 除 以 一 个 1/ VN 得 到 传 里 叶 基 F: 
本 文 利用 相关 性 能 优异 的 标准 正 交 傅 里 叶 基 , 构造 P-Ifourier 矩 ri 1 1 1 
阵 作为 观测 矩阵 。P-Ifourier 矩阵 的 构造 方法 如 下 : 1 W Ww We” 
a) 产生 一 个 传 里 叶 变 换 和 矩阵 Fx，P,eC™*，Cw*w 是 N r= Wi wt Wa (10) 
XN 维 复数 域 空间 。 | We Wa Wa 
1 1 1 1 1 
1 Ww WwW? WwW! 于 
Fs=1l Wi Ws Wv (9) 
] Wi Wan | Wa Den) 
FF = 下 
1 1 1 1 1 1 1 EE 1 
D> Be 0 i a Ca 
1 ew ew e 1 ew ew ee 
—(N-D)27j i -CXN-DCN-D.2z7 CX-D.2z7 ee (CN-DCN-D22z/ 
11 e » e 人 .ee a [ee BA e 人 | 
7 为 NXN 维 单位 矩阵 ， 因 此， 和 矩阵 是 酉 矩阵， 在 西 空 
间 内 是 标准 正 交 浊 5 之 (14) 
c) 对 傅 里 叶 基 瓦 求 着 得 到 了 逆 傅 里 叶 基 (D: 
上 为 半径 的 圆 盘 。 可 知 
| (mi) (wa) (wy 4 CA A) hm A) =1 (15) 
| ， 根据 式 (7)(15) 可 以 得 出 , P-Ifourier 矩阵 4 不 仅 满足 RIP 特 
7 rl Wl (we™™) 性 ， 而 且 其 RIP 特性 比较 理想 。 
1 (wi™) (We )) . (Wwe) 3 ”仿真 分 析 
be es i 为 了 验证 本 文 所 提出 的 频谱 感知 模型 的 可 行 性 ， 以 及 所 构 
I de Kee 建 的 P-Ifourier 观测 矩阵 A 的 重 构 性 能 以 及 重 构 精度 相 比 于 高 
= 和 斯 随机 观测 矩阵 和 文献 3] 所 提出 的 确定 性 嵌入 式 混沌 序列 - 循 
: : 环 Toeplitz 结构 的 观测 矩阵 是 否 有 显著 改善 ， 本 节 利 用 
ee MATLAB 仿真 软件 进行 实验 分 析 。 
: 3.1 ”模型 可 行 性 验证 实验 
d) 对 逆 傅 里 叶 基 g 随 机 选取 M 行 向 量 便 构建 了 P-Ifourier 本 文选 用 的 时 域 干扰 信号 x 为 多 音 干扰 ， 采 样 频率 512， 
矩阵 4。P-Ifourier 矩阵 4 保持 了 傅 里 叶 基 的 标准 正 交 特性 ， ”干扰 信号 x 的 频 域 稀 琉 度 K=7, 稀疏 基 顾 为 标准 正 交 侍 里 叶 基 ， 
并 且 jx(A) =1， 因 此 具有 理想 的 非 相 关 性 ， 下 面 就 4 是 否 满足 。 观测 矩阵 为 本 文 所 构建 的 P-Ifourier 观测 矩阵 ,观测 维 数 M=128， 
RIP 特性 进行 证 明 。 恢复 算法 选择 OMP 算法 ， 仿 真 结果 如 图 2~4 所 示 。 
2.3 ”观测 矩阵 的 RIP 特性 证 明 图 2 为 待 检 测 的 时 域 多 音 干扰 信号 x, 可 以 看 出 , x 在 时 域 
设 矩 阵 大 = 4a4 ,由 2.2 节 可 知 上 并 不 具有 稀 琉 特性 ， 因 此 利用 标准 正 交 傅 里 叶 基 屎 对 x 进行 
De (13) 各 邪 表达 ， 图 3 表示 的 是 PIfourier 观测 算 阵 根据 * 的 频 域 和 
为 =0 (iz#) 玻 特性 而 进行 压缩 采样 后 的 采样 波形 。 根 据 图 2 和 图 3 可 以 看 
根据 Gersgorin 国 盘 定理 , 即 <p ，D; 是 复 平面 上 以 霹 。 出 ”采样 过 程 并 非 是 依据 信号 的 带宽 而 进行 的 均匀 采样 ， 而 是 
站 根据 干扰 信号 频 域 的 结构 特征 而 进行 的 非 均匀 采样 ， 使 信号 的 
为 圆心 ， 以 采样 频率 降低 了 75%。 
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图 4 为 x 的 原始 频谱 与 基于 P-Ifourier 观测 矩阵 恢复 的 频 


谱 对 比 示 意图 ， 


此 可 以 直观 地 表明 本 文 所 构建 的 P-Ifourier 观 


测 和 矩阵 能 够 够 较为 精准 的 恢复 
本 文 所 提出 模型 的 可 行 性 。 


出 x 的 频谱 。 一 定 程度 上 验证 了 


dp 
-| | | | i | | 


0.5 
时 间 /s 
图 2 时 域 多 音 干扰 
80 
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图 3 信号 压缩 采样 波形 
25F- 
2 一 恢复 频谱 
立 原始 信号 的 频谱 
20 上 上 华 


图 4 


3.2 重 构 性 能 分 析 


上 一 小 节 从 比较 直观 的 角度 验证 了 模型 的 可 行 性 ， 本 小 节 
将 对 P-Ifourier 观测 矩阵 的 重 构 性 能 与 高 斯 随机 观测 矩阵 和 文 
献 [3] 所 提出 的 确定 性 峰 


恢复 的 频谱 与 信号 原始 频谱 的 对 比 


入 式 混 沌 序列 -循环 Toeplitz 结构 的 观 


测 矩 阵 进行 对 比分 析 ， 待 检测 信号 带宽 N=512， 稀 玻 度 玉 的 变 


化 范围 为 S~70， 不 考虑 噪声 影响 ， 观 测 维 数 M=128， 重 构 算法 


选取 OMP 算法 ， 


蒙特 卡 洛 仿真 10000 次 。 仿 真 结 果 如 图 $ 所 


刘 洋 ， 


不 。 


: 基于 P- jo hin 算 阵 的 完 


inaXiv 合 作 期 刊 


压缩 感 知 i 


重 构 概率 /% 


图 5 无 噪 条 件 下 观测 矩阵 重 构 概率 对 比 


5 为 无 噪 条 件 下 不 同 观 测算 阵 的 重 构 概率 对 比 示 意图 。 


图 可 知 ， 在 无 噪 条 件 下 ， 当 信 
观测 矩阵 都 能 以 接近 100% 的 概率 台 


列 -循环 Toeplitz 
Ifourier 观测 和 矩阵 仍 能 


列 -循环 Toeplitz 


在 25~45 范围 内 时 ， 高 斯 随机 观测 矩 
结构 的 观测 矩阵 的 习 
以 接近 100% 的 概率 重 构 原始 信号 


号 稀 疏 度 玉 小 于 25 时 ， 三 种 
构 原 始 信号 


; 当 稀 疏 度 天 
阵 和 确定 性 嵌入 式 混沌 序 
E 构 概率 急剧 下 降 ， 而 P- 
; 当 稀 


疏 度 K 高 于 50 时 ， 高 斯 随机 观测 算 
结构 的 观测 矩阵 已 无 法 重 构 原始 信号 ,信号 失 


阵 和 确定 性 嵌入 式 混沌 序 


号 ,之 后 才 缓慢 下 降 。 
然而 现实 生活 中 的 感知 环境 通常 


真 , 而 P-Ifourier 观测 矩阵 在 K=50 时 仍 能 以 98% 的 概率 重 构 信 


是 较 低 信 噪 比 的 环境 ， 接 


下 来 进行 加 噪 实验 , 信 品 比 为 5dB, 稀 疏 度 K 变 化 范围 为 5~45， 


其 余 条 件 同 上 ， 念 真 结果 如 图 6 所 示 。 


~ AR 


中 SS 、 


* 一 高 斯 随机 观测 矩阵 
嵌入 式 混沌 序列 -循环 Toeplitz 结 构 的 观测 矩阵 
专 一 P-IFourier 观 测 矩 阵 
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图 6 ， 信 噪 比 为 dB 条 件 下 观测 矩阵 重 构 概 率 对 比 


6 表明 在 较 低 信 噪 比 条 件 下 观测 矩阵 的 重 构 性 能 : 
内 时 ， 三 种 矩阵 都 入 


定 程 度 的 下 降 。 当 天 在 5~8 的 范围 


以 接近 


100% 的 概率 完成 信号 重 构 ; 当 
观测 拢 


循环 Toeplitz 


阵 和 确定 性 嵌入 式 混 沌 序列 -循环 Toeplitz 
阵 的 重 构 性 能 急剧 下 降 ， 在 K=20 时 ， 确 定性 嵌入 式 混 沌 序列 - 
结构 的 观测 矩阵 的 重 构 概 率 为 19.5%， 高 斯 随机 


天 在 8~20 范围 内 时 ， 高 斯 随机 


结构 的 观测 甜 


观测 矩阵 的 重 构 概 率 为 17.1%， 无 法 满足 重 构 要 求 ,信号 基本 失 


真 ， 而 P-Ifourier 观测 矩阵 仍 能 


以 96.1% 的 概率 重 构 信 号 ; 当天 
大 于 20 时 了 P-Ifourier 观测 矩阵 的 重 构 性 能 才 缓慢 下 降 。 
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3.3” 重 构 精 度 分 析 
本 次 实验 的 信 噪 比 变化 范围 为 -10~-1dB, 其 余 条 件 与 3.1 节 


相同 , 重 构 精度 的 分 析 指 标 为 均 方 误差 (mean square error, MSE): 


Es 
MSE=— 2 (16) 
E | 
下 
仿真 结果 如 图 7 所 示 。 
J 
一 二 一 P-IFourier 观 ; FF 
i -循环 Toeplitz 结 构 观测 矩阵 
矩阵 
本 
ss 
忌 
10o4 Re 
， 
E k 
1002 上 
10 9 8 7 | 5 4 3 


信 品 比 /dB 


图 7 不 同 信 噪 比 下 观测 矩阵 的 均 方 误差 


| 


给 出 了 各 观测 矩阵 在 不 同 信 品 比 下 的 均 方 误差 ， 由 图 
可 知 ， 各 个 观测 矩阵 的 均 方 误差 随 信 噪 比 的 增加 而 降低 ，P- 
Ifourier 观测 珑 阵 的 均 方 误差 曲线 明显 优 于 高 斯 随机 观测 珑 阵 
和 确定 性 岩 入 式 混沌 序列 -循环 Toeplitz 结构 的 观测 矩阵 的 均 方 
误差 曲线 , 这 是 由 于 P-Ifourier 观测 矩 阵 的 原始 方 阵 在 西 空间 内 
是 标准 正 交 和 矩阵 ， 因 此 具有 理想 的 非 相关 特性 ， 而 其 他 两 种 观 
测 矩 阵 的 相关 性 与 前 者 相 比 存在 一 定 的 不 足 。 经 计算 ,P-Ifourier 
观测 矩阵 的 均 方 误差 比 确定 性 嵌入 式 混 沌 序列 -循环 Toeplitz 结 
构 的 观测 矩阵 平均 降低 3.1%, 比 高 斯 随机 观测 矩阵 平均 降低 约 
4.7%。 因 此 本 文 构造 的 P-Ifourier 观测 矩阵 在 重 构 精度 方面 相 
比 于 其 他 观测 矩阵 有 较为 明显 的 改善 。 


| 


TOI 


4 ”结束 语 


压缩 感知 理论 应 用 于 传统 的 频 六 
号 的 采样 频率 ， 但 是 在 重 构 性 能 和 重 构 精度 方面 存在 一 定 的 不 
足 。 本 文 利用 正 交 性 能 和 非 相 关 性 能 理想 的 标准 正 交 健 里 叶 基 
构造 观测 矩阵 ， 并 将 该 矩阵 应 用 于 宽带 压缩 频谱 感知 模型 中 ， 
该 矩阵 能 保持 标准 正 交 传 里 叶 基 的 相关 性 能 ， 在 重 构 性 能 和 重 
构 精 度 方面 明显 提高 ， 是 一 种 性 能 优良 的 频谱 感知 观测 矩阵 。 


普 感知 问题 能 够 有 效 降 低 信 


hinaxiy 合 作 期 于 


刘 洋 ， 等 : 基于 了 -Ifourier 和 纶 阵 的 宽带 压缩 有 
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